Static impurities in the kagome lattice: dimer freezing and mutual repulsion 
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We consider the effects of doping the 5 = 1/2 kagome lattice with static impurities. We demon- 
strate that impurities lower the number of low-lying singlet states, induce dimer-dimer correlations 
of considerable spatial extent, and do not generate free spin degrees of freedom. Most importantly, 
they experience a highly unconventional mutual repulsion as a direct consequence of the strong 
spin frustration. These properties are illustrated by exact diagonalization, and reproduced to semi- 
quantitative accuracy within a dimer resonating-valence-bond description which affords access to 
longer length scales. We calculate the local magnetization induced by doped impurities, and consider 
its implications for nuclear magnetic resonance measurements on known kagome systems. 

PACS numbers: 75.10.Jm, 75.30.Hx, 76.60.-k 



I. INTRODUCTION 

The kagome geometry (Fig. 1) presents one of the most 
highly frustrated quantum spin systems achievable in 
two dimensions (2d) with only nearest-neighbor Heisen- 
berg interactions. The S — 1/2 kagome lattice has been 
found to have a spin-liquid ground stated whose ultra- 
short spin-spin correlation lengths are effectively those 
of a dimer liquid; 2 ^ and whose excitation spectrum- 
shows a manifold containing an extremely large num- 
ber of low-lying singlets^ These properties are repro- 
duced very well by a short-range resonating-valence- 
bond (RVB) description^ based in fact only on nearest- 
neighbor dimer formation (Fig. 1). 

On the theoretical level, a similar degree of frustration 
is found only in quantum dimer models, 7,8 which have 
also been considered recently in the kagome geometry? 9 , 
and in models containing higher-order spin terms and 
multiple-spin exchange. 1 - Experimentally, a variety of 
materials displaying the kagome structure is known to 
exist, although to date none have been found which con- 
tain an ideal, 2d, vacancy-free system of spins S = 1/2. 
These include the jarosites (H 3 0)Fe 3 (OH) 6 (S0 4 )2ri 
KFe 3 (OH) 6 (S0 4 )2, and KCr 3 (OH) 6 (S0 4 )2±2*±2 and the 
volborthite Cu 3 V 2 7 (OH)2.H 2 Oi4 The best charac- 
terized of these compounds is the magnetoplombitc 
SrCrg p Gai2_9 P Oi9 (SCGO), 15 in which some of the S 
= 3/2 Cr 3+ ions form approximately 2d kagome planes. 
However, because in this case the ideal stoichiometry 
(p = 0) remains unachievable, discussion of the influ- 
ence of static (and also, for Ga 3+ , spinless) impurities 
represents an important issue. 

Studies of static vacancies in spin liquids have tended 
to focus on the question of induced free spin degrees of 
freedom— In systems such as spin ladders, and also for 
the 2d square lattice^ 7, the total energy appears to de- 
crease as such holes in the spin background are made to 
approach each other, implying an effective mutual attrac- 




FIG. 1: The kagome lattice. Thick black lines represent an ar- 
bitrarily chosen covering by nearest-neighbor dimers formed 
between S — 1/2 spins on each lattice site. 



tion. It is generally believed that in an RVB description 
of the ground state the minimal disruption of the wave 
function, and thus the most favorable energy state, oc- 
curs when holes occupy adjacent sites, thus generating a 
simple pairing mechanism for mobile dopants. The issue 
of impurites in frustrated magnets has received less atten- 
tion, and appears to be more involved: it has been pro- 
posed that an absence of induced free spins corresponds 
to deconfined, spinon-like excitations: 18 Consideration of 
the problem for the Id frustrated chain; 1 " 9 which provides 
an example with no free spins localized around the va- 
cancies, suggests the importance of degenerate singlets 
in the ground-state manifold. This result indicates that 
the kagome geometry presents the most probable can- 
didate for exotic behavior— as a consequence of strong 
frustration in 2d. 

In this study we investigate the effects of static impu- 
rities in the kagome lattice. We find that, despite the 
very short intrinsic spin-spin correlation lengths, dimer- 
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dimer correlations develop over considerable distances in 
the presence of dopants. Further, we show that there is 
no evidence for free local moments induced around im- 
purity sites, and that a highly unusual effective repulsive 
interaction arises between such holes in the spin back- 
ground. 

In Sec. II we provide an outline of the numerical and 
anaytical methods to be employed, and discuss the fun- 
damental effects on the excitation spectrum of doping by 
spinless impurities. In Sec. Ill we consider the case of one 
hole in an odd-site cluster, demonstrating "dimer freez- 
ing" in the spin-spin correlation functions and dimer- 
dimer correlations of anomalously long range. The in- 
duced magnetizations in the spin sectors S — 1 and higher 
are used to illustrate the consequences for nuclear mag- 
netic resonance (NMR) measurements. Sec. IV presents 
spin-spin correlation functions which characterize the 
situation with two impurities, supplemented by total- 
energy analyses quantifying the novel repulsive interac- 
tion between the holes which emerges from this study. 
Sec. V contains further discussion and our conclusions. 
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FIG. 2: Integrated density of states for an 18-site S — 1/2 
kagome cluster, computed by exact diagonalization both with- 
out impurities (solid line) and with 2 impurities (dashed line). 
The vertical lines mark the position of the first triplet state 
in the absence of impurities (solid) and in their presence 
(dashed). 



II. ELIMINATION OF SINGLET STATES AND 
ABSENCE OF LOCALIZED SPINS 

A. Exact Diagonalization 

By "impurity" we refer henceforth to nonmagnetic im- 
purities. The techniques we employ are based primarily 
on numerical calculations for small clusters with periodic 
boundary conditions chosen as symmetrically as possi- 
ble. Specifically, we perform exact diagonalization (ED) 
studies of the Heisenberg Hamiltonian 

F = J^S r S,, (1) 

m 

where J is the antiferromagnetic superexchange interac- 
tion and (ij) denotes nearest-neighbor sites, for clusters 
up to 27 sites with one impurity, and up to 24 sites with 
two impurities for all available configurations. We ana- 
lyze the total energies, the energies per doped impurity, 
and also the form of the energy spectra, focusing on the 
number of singlet states in the low-energy manifold be- 
low the first triplet for the finite systems under consider- 
ation, as well as on the number of these states removed 
by impurity doping. The frustrated spin interactions in 
the presence of impurities are characterized for the entire 
cluster by the bond spin-spin correlation functions and 
by the induced site magnetizations in every spin sector. 

The density of states integrated over the lowest-lying 
energy levels is shown in Fig. 2 for systems with and 
without impurities. Fig. 2 compares ED results from 
systems of 18 sites, for which in the impurity configura- 
tion illustrated 8 of the 12 singlet states occurring below 
the first triplet, indicated by the vertical lines, are elim- 
inated by the presence of 2 impurities. Although the 



separation of the triplet from the lowest remaining sin- 
glet decreases, this energy gap always remains a finite 
fraction of the superexchange interaction J. Thus the 
most important qualitative result emerging even for the 
smallest systems is that, in contrast to spin-ladder and 
dimerized or quadrumerized systems, no low-lying triplet 
state is induced by the presence of impurities. This im- 
plies directly 19 that there are no quasi-free spin degrees 
of freedom, by which is meant spins effectively isolated 
from their local environment and interacting with each 
other only on an effective energy scale much smaller than 
J. The fact that no localized moments are formed in the 
vicinity of impurities may be ascribed to the continuum 
of singlets beyond the ground state, which has the same 
function as the two-fold degeneracy of the ground state 
in the frustrated chain ^ namely that the rearrange- 
ment of dimers incurs no energy penalty. The absence 
of local spins around nonmagnetic impurities represents 
one of the most significant differences between frustrated 
and unfrustrated systems, and has been further inter- 
preted as implying a deconfinement of the spin degrees 
of freedom^ 

Table 1 provides more details of the number of states 
in the singlet manifold which are removed by impurity 
doping for clusters of 12 and 18 sites, for all possible 
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12 2 3 
4 4 4 4 
2 2 3 1 


1 2 2 3 3 
12 12 12 12 12 
3 4 4 3 4 



TABLE I: Number of low-lying singlets below the first triplet 
for kagome clusters without impurities (N°, third row) and 
with two impurities (iVf, fourth row) as a function of the 
Manhattan distance <1m between impurities. Results obtained 
by ED are shown for clusters of N = 12 and 18 sites. 
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FIG. 3: Integrated density of singlet states for an 18-site S = 
1/2 kagome cluster, calculated in the dimer-RVB basis both 
without impurities (solid line) and with 2 impurities (dashed) . 



configurations of 2 impurities. The numbers vary little 
as a function of the proximity of the impurities, and the 
average number of low-lying states which are eliminated 
is in agreement with simple considerations of the num- 
ber of singlet states in the small system (below) . We have 
also performed analogous calculations for 24-site clusters 
of different shapes. For all clusters the same behavior 
is found upon introducing two impurities, in that the 
number of singlets below the first triplet is reduced, but 
here the lowest-lying excitation is not a triplet. However, 
these clusters do not appear to be representative of the 
generic kagome physics, as the number of low-lying sin- 
glets in the absence of impurities is significantly smaller 
than the expected valued of 1.15 N . This situation is pre- 
sumed to be a consequence of the low symmetry of the 
24-site clusters, some of which are essentially Id systems, 
and as a result the densities of states are not presented 
here. 



B. Dimer-RVB Basis 

We calculate in addition the ground-state properties, 
namely the energy and spin-spin correlation functions, of 
clusters doped with one and with two impurities within 
the dimer-RVB basis. This is the basis of states in which 
every spin forms a singlet dimer with one of its nearest 
neighbors in such a way that every spin of the lattice is 
involved in a dimer (Fig. 1). The use of this basis& natu- 
rally implies a truncation of the total number of available 
states in the manifold of low-lying singlets, but we will 
show below that it provides a satisfactory description of 
the ground state of the doped system. Dimer-RVB cal- 
culations permit the consideration of somewhat larger 
cluster sizes, and we will present results up to and in- 
cluding 39 sites with one impurity and 36 sites with two 
impurities. 

Calculations performed in the dimer-RVB basis yield 



results very similar to the exact ones for the integrated 
density of states (Fig. 3), giving clear evidence of a strong 
linear correlation (see also Fig. 9) between the ED and 
dimer-RVB spectra in both 0- and 2-impurity cases. 
Demonstration of a quantitative correspondence between 
the spectra would require shift and scaling factors aris- 
ing from an energetic renormalization related to the fact 
that calculations in the dimer-RVB basis are variational 
in nature, so that the true eigenstates are dressed RVB 
wave functions, but this procedure is not necessary for 
the analysis to be presented here. 

Given that the RVB basis represents a very signifi- 
cant truncation of the total Hilbert space, the qualita- 
tive similarity of the dimer-RVB spectra to those ob- 
tained from the exact calculations, both with and with- 
out impurities, provides important information about the 
nature of the exact wave functions of the kagome sys- 
tem. This observation is further confirmed by the semi- 
quantitative agreement with exact results for spin cor- 
relation functions and ground-state energies computed 
using the dimer-RVB basis. We defer this comparison to 
Sees. Ill and IV, and comment here only that our results 
in the presence of impurities provide additional strong 
justification for the statement that the ground state of 
the kagome spin system is described appropriately by a 
maximally short-ranged RVB wave function, presumably 
as a consequence of the very strong frustration. 

The dimer-RVB basis is a natural framework for un- 
derstanding the absence of localized spins around an im- 
purity site, as all remaining spins are simply incorporated 
within rearranged singlet dimer coverings. It also pro- 
vides a qualitative reflection of the "evaporation" of low- 
lying singlets from the ground-state manifold (Figs. 2-3 
and Table I) which arises on doping a cluster with two 
impurities, and in fact yields additional insight into the 
mechanism for this exclusion. It can be shown that the 
total number of singlet dimer coverings of a doped system 
is given by 



coverings 



_ 2N/3-N imp +N G + N c 



(2) 



where N is the number of sites of the system under con- 
sideration, N- lmp the number of doped impurities, Nq 
the number of independent clusters into which the sys- 
tem is divided by the dopant distribution, and iV e t the 
number of triangles containing 3 impurities. A complete 
derivation of this result lies beyond the scope of the cur- 
rent analysis, and we present only the following heuristic 
explanation. By reducing each triangle to an effective 
S = 1/2 degree of freedom;^ the dimer basis involves a 
truncation of the 2 N total states to 2 Ar / 3 , multiplied by 
an additional factor of 2 for the overall spin direction on 
each independent cluster. Each added impurity removes 
one spin, halving the number of states in every mani- 
fold and thus justifying the dependence on N- lmp . The 
only exception to this statement occurs when the added 
impurity falls on a triangle which already contains two 
impurities, as a result of which it is already fully con- 
strained and removal of the third site no longer affects 
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FIG. 4: Spin-spin correlations for a single impurity in a 
27-site cluster, obtained by ED. Bar width represents the 
strength of correlation functions on each bond on a linear scale 
where the strongest correlation function is (Sj-Sj) = —0.69. 
Dashed lines denote bonds on which (Sj-Sj) > 0. 



the number of coverings. Eq. can be shown to remain 
valid for all vacancy concentrations, including most no- 
tably values of JVi mp exceeding the percolation threshold, 
beyond which Nq > 1, under the condition that the im- 
purity distribution does not cut the system into segments 
containing odd numbers of spins. 



III. SPIN CORRELATION FUNCTIONS AND 
INDUCED MAGNETIZATION 

A. Dimer "freezing" 

In this section we consider the effects of a single impu- 
rity introduced in an odd-site cluster. The bond spin- 
spin correlation functions (Sj-Sj) around a single impu- 
rity may be computed by ED for clusters of up to 27 
sites. The ED calculations permit in addition the de- 
termination of the magnetizations induced on each site 
by the presence of the impurity, with results for every 
spin sector corresponding to the situation in all applied 
magnetic fields, which are presented in Subsec. C below. 

The correlation functions obtained by ED are illus- 
trated in Fig. 4 for a 27-site cluster with periodic bound- 
ary conditions. The most striking feature of these results 
is that dimer bonds are effectively frozen in the two de- 
pleted triangles, taking a value (S, • Sj) ~ —0.69 close 
to the theoretical minimum of —3/4 for a pure singlet 
state. For comparison, the average value of (Sj • Sj) in 
the pure system, which retains full freedom to resonate, 
is a factor of 3 smaller (—0.22)^ The bonds neighboring 
the frozen dimers then show very weak correlations, while 
the next neighbors are again enhanced, and the spin cor- 
relations oscillate with distance away from the impurity 
[also illustrated in Fig. 5(b)]. Thus although spin-spin 
correlations in the kagome system are extremely short- 
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FIG. 5: (a) Spin-spin correlation functions G(cLm) = (Si ■ Sj) 
on nearest-neighbor bonds as a function of their Manhattan 
distance <1m from a single impurity, for clusters of 21 and 27 
sites, (b) Spin-spin correlations as a function of (Manhattan) 
bond distance from a single impurity, obtained by dimer-RVB 
calculations for a 39-site cluster. Points represent the average 
over all bonds at the same distance from the impurity, the 
solid line the value (Sj-Sj) = —0.22 for the impurity-free 
system. 



ranged, 2 induced dimer-dimer correlations show a con- 
siderable spatial extent, reaching in fact well beyond the 
sizes obtainable by studying small clusters using either 
ED or dimer-RVB techniques. For the limited data avail- 
able it is unfortunately not possible to specify the func- 
tional form of the decay to equilibrium. Further analysis 
of this issue, based on data obtained from mutual hole 
repulsion, is presented in Sec. IV. 



B. Dimer-RVB Basis 

Fig. 5(a) compares the bond spin correlation values 
(Sj • Sj) obtained from ED and from the dimer-RVB ba- 
sis for all bonds in clusters of 21 and 27 sites with one 
impurity. These clusters represent the maximal sizes ac- 
cessible by ED, and for such systems it is clear that the 
dimer-RVB basis reproduces the exact values with quan- 
titative accuracy. From the nature of the RVB basis it 
is perhaps to be expected that nearest-neighbor bond 
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FIG. 6: Site magnetization profile of the lowest triplet state 
induced by a single impurity in a 27-site cluster, obtained 
by ED. Dot radius represents the magnitude of the local mo- 
ment on a linear scale where the largest dot corresponds to 
a moment of 0.13/xb, and the grey dots represent sites with 
induced moments opposite to the field. 



correlations are treated more accurately than any other 
physical property. We make use of this feature to extend 
the calculations of spin correlation functions to a system 
of 39 sites with oie impurity. Because on this cluster 
the boundary conditions lead to a rather low symmetry, 
Fig. 5(b) shows the value of (Sj-Sj) averaged over all 
bonds at the same Manhattan distance from the impu- 
rity, which illustrates both the oscillatory nature of the 
frozen dimerization pattern and the very long range over 
which an impurity acts to disrupt the spin configuration. 



C. Induced magnetic moments 

The local magnetization pattern induced in the vicin- 
ity of a doped nonmagnetic impurity may be computed 
by exact diagonalization in the sector of total spin (S 1 , S z ) 
— (1,1), and is displayed for the 27-site cluster in Fig. 6. 
The site magnetizations exhibit a certain anticorrelation 
with the strength of the dimers formed on each bond, and 
thus have a tendency to oscillate with increasing distance 
from the impurity. The induced moments are particularly 
enhanced on sites far from the impurity, where different 
paths on the periodic cluster interfere, and the interfer- 
ence in this regime may even result in moments oriented 
opposite to the effective magnetic field (Figs. 6, 7). The 
magnetization profile of the surrounding sites is the sin- 
gle most important consequence of impurities which is 
reflected in local-probe measurements of doped spin sys- 
tems. 

Fig. 7 shows the magnetization per site induced by 
proximity to a nonmagnetic impurity for each spin sec- 



FIG. 7: Induced magnetization per site {ii as a function of 
the spin sector S z (analogous to an applied magnetic field) 
for all inequivalent sites— neighboring an impurity, computed 
by ED for a 27-site cluster. With reference to the site labels 
in Fig. 6, the separate lines i represent 1st (4, 5, 23, 24), 2nd 
(6, 9, 22, 25), 3rd (10, 12, 19, 20), 4th (2, 3), equal 5th [(7, 
8, 26, 27) and (13, 15, 16, 18)], 7th (14, 17), and 8th (11, 21) 
neighbors of the impurity site (1). 



tor obtainable by ED for a 27 site cluster. The spin 
sector may be considered as equivalent to an effective 
magnetic field which varies between zero (S z = 0) and 
saturation (S z = 13). The most pronounced variation 
is observed for the nearest-neighbor sites (filled circles in 
Fig. 7), which have very small magnetization at low fields 
because they are bound into the strong dimers frozen by 
the presence of the impurity on the same triangle. How- 
ever, between S z = 4 and S z — 6 the moment shows an 
abrupt jump to saturation, implying the breaking and 
full polarization of these dimers at fields in excess of 1 /3 
of the saturation field. We speculate that this behavior 
may be due to the presence of a particularly stable state 
of the remaining spins which corresponds to a 1/3 mag- 
netization plateau of the undoped kagome latticed The 
moments of the further-neighbor sites show a generally 
less systematic evolution, usually with a non-monotonic 
character which implies a very significant rearrangement 
of the induced local magnetization as a function of field. 
This evolution is again marked by the most significant 
changes occurring close to the net magnetization value 
of 1/3. 

NMR has emerged as the technique of choice for an- 
alyzing in real space the effects of doped impurities on 
magnetizations and spin correlations. Information con- 
cerning these quantities is extracted from Knight-shifts, 
linewidth alterations and, for certain cases, from re- 
laxation times. A detailed NMR and NQR (nuclear 
quadrupole resonance) analysis has been conducted re- 
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FIG. 8: Energy A2 for a pair of holes as a function of hole 
separation, obtained by ED for clusters of 12, 18 and 24 sites 
with two impurities. 

cently for 71 Ga impurities in the kagome planes of the S 
= 3/2 spin system SrCr 9 pGa 12 _ 9p Oi9 (SCGO) While 
a quantitative determination of any shifts or linewidths 
depends on a detailed Hamiltonian for the hyperfine in- 
teractions in a specific material, the results obtained in 
Fig. 7 may be considered as illustrative for the kagome 
geometry. 

The important qualitative features are the absence of a 
single spin-1/2 degree of freedom anywhere in the vicin- 
ity of the impurity, and the rather large extent of the area 
of affected sites in spite of the very short spin correlation 
lengths. The latter effect is contained in Fig. 7 in that 
the induced site magnetization values on the vertical axis 
are directly proportional to the peak shifts measured in 
an NMR field sweep. Impurities in the kagome lattice 
therefore produce a large number of satellite peaks. 

Finally, with regard to the fact that most known 
kagome systems, including SCGO, are not S = 1/2, there 
remain rather few studies of this geometry for higher val- 
ues of the spini22iS On rather general grounds based on 
the geometry of the lattice, one expects that both tri- 
angles and depleted triangles retain the integral or half- 
integral nature of the spin per site. Thus it is at least 
plausible that the results for the 5 = 1/2 system are di- 
rectly relevant to higher half-integral-spin systems such 
asSCGO^ 



IV. HOLE REPULSION 

In this section we present results characterizing the be- 
havior of an even-site cluster doped with two impurities. 
A comparison with the undisturbed system^ reveals that 
the average energy per remaining site not only increases 
in the presence of impurities, but increases further as 
the impurities are brought closer together. Both of these 
features are evident from the ED results of Fig. 8, which 
shows the energies of all different two-hole states avail- 
able for the relevant system sizes. The higher energy 
per site of doped systems confirms, by considering the 



FIG. 9: Comparison of energies A2,ed and A2,rvb obtained 
respectively from ED and from calculations in the dimer-RVB 
basis for a pairs of holes of different separations on clusters of 
12, 18, and 24 sites with two impurities. 

contributing processes, that the net loss of resonance en- 
ergy for dimers constrained by the presence of impurities 
to occupy only the bonds on which they are essentially 
"frozen" (Sec. Ill) is not compensated by the gains in 
dimerization energy for those bonds which are enhanced. 

The effective hole repulsion is most unusual in low- 
dimensional spin systems ^L2& and in fact we are unaware 
of analogous results for any model other than the frus- 
trated chaini^i This repulsion of holes constitutes the 
central result of our analysis: particularly in the con- 
text of short-range RVB wave functions, it contradicts 
directly the most straightforward expectation for the be- 
havior of doped holes, namely that when mobile these 
would choose to minimize the net disruption of the sys- 
tem by occupying the same dimer, and thus would ex- 
hibit a pairing attraction^ The consequences for mobile 
impurites with repulsive interactions may be expected 
to include hole crystallization, phase separation and the 
possibility of transitions between metallic and insulating 
phases, and would bear considerable further investiga- 
tion. 



A. Dimer-RVB Basis 

We illustrate in Fig. 9 the efficacy of the dimer-RVB 
basis for the purposes of comparing ground-state ener- 
gies. The correspondence between the dimer-RVB and 
ED results is shown by comparing the energies of all two- 
hole states in clusters of 12, 18 and 24 sites with two 
impurities, which represent all system sizes accessible by 
ED. Full details of the extraction of the energy values A2 
in Fig. 9 are given below. Although an inspection of the 
abcissae reveals that the agreement between individual 
pairs of values is not as convincing as for the spin cor- 
relation functions (Sec. Ill), the main difference between 
the two sets of data may be reduced to a relative upward 
shift of all the RVB energies. The direct linear correla- 
tion between ED and dimer-RVB results (alluded to in 
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FIG. 10: Spin-spin correlations for two impurities in a 36- 
site cluster, obtained in the dimer-RVB basis. Bar width 
represents the strength of correlation functions for each bond 
on a linear scale where the largest value is —0.66. 
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FIG. 11: Energies A2 for two-hole states as a function of hole 
separation for all configurations of two impurities on 30- and 
36-site clusters, obtained by dimer-RVB calculations. The 
solid line marks twice the single-hole energy 2Ai relevant for 
comparison, taken from the 33-site cluster. 



Sec. II) demonstrates unambiguously that comparing the 
energies of different configurations of the two impurities, 
which is the primary goal of the energy analysis to follow, 
is a valid qualitative exercise when performed using the 
dimer RVB results. 

We show first the spin-spin correlation functions in the 
presence of a pair of impurities (Fig. 10). When consider- 
ing the overlap of single-impurity dimerization patterns 
(Sec. Ill), it is evident that these are not fully mutu- 
ally compatible. In fact, one finds by inspection that not 
only the configuration illustrated in Fig. 10 but all pos- 
sible configurations of two impurities on the kagome lat- 
tice are in some measure incompatible. Comparison with 
other quantum spin systems suggests that this complete 
incompatibility is very specific to the kagome geometry, 
and its extent may be quantified by extracting the effec- 
tive interaction energy of the holes (Fig. 11). In cluster 
calculations the frustration is also shown by the presence 
of relatively large local moments induced on the sites at 
the mismatch between dimerization regimes (cf. Fig. 6). 



B. Interaction between impurities 

Fig. 11 shows twice the energy per hole as a function of 
separation for all configurations of two impurities on 30- 
and 36-site clusters, obtained in the dimer-RVB basis. 
We have analyzed the data in the following manner. E n 
is the ground-state energy of each cluster in the absence 
of impurities, and is computed for even clusters then ex- 
trapolated for odd clusters from the results for systems 
of sizes N — 3 and N + 3. E\ is the ground-state energy 
in the presence of a single impurity, which is computed 
for odd clusters and may then be extrapolated for even 
clusters from the results for systems of sizes N — 3 and 
iV+3. Ei is the ground-state energy with two impurities, 
and again may be computed directly for even clusters and 



extrapolated for odd clusters. We then define 

Ai = Eh. - E (3) 

as the energy of one hole, 

A 2 = E 2 - E (4) 

as the energy of a two-hole state, and 

e = A 2 - 2Ai = E 2 - 2Bh + E (5) 

as the energy of interaction between two holes. Because 
of the extrapolation procedure it is important that the 
data sets for different system sizes be comparable. While 
certain problems arise for some system sizes due to the 
relatively large influence of the boundary conditions, we 
find the process to be valid for systems in excess of 27 
sites. In this connection we note that the two-hole ener- 
gies A 2 shown in Fig. 11 are independent of the extrap- 
olation between different system sizes, and it is only the 
direct comparison with Ai for the 33-site system which 
requires validation of the procedure. 

It is clear that the effective mutual repulsion decays 
only slowly. Because of the restricted sizes of the acc- 
cessible clusters we are unable to speculate on the func- 
tional form of this decay, and in particular on the pos- 
sibility that it may be long-ranged. We note only that 
for short distances (3-4 lattice constants) it appears to 
exhibit an approximate \jd r dependence, beyond which 
a steeper decrease is observed which may be either an 
intrinsic property of the interaction or a consequence of 
the system size. The form of the decay appears similar to 
that of the bond spin correlation functions towards their 
undisturbed value [Fig. 5(b)]. 

At the phenomcnological level, a repulsive interaction 
between two impurities may be regarded as surprising 
in that their interference might have been expected to 
suppress the freezing of individual dimer patterns, lead- 
ing to a restoration of resonance energy. Our findings to 
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the contrary indicate that the dominant effect of increas- 
ing increasing impurity separation lies in the reduction in 
strength of the partially frozen single-impurity dimeriza- 
tion patterns at the boundaries where they are forced to 
mesh, which reduces the corresponding energetic penalty. 
Given the slow decay of the dimer-dimer correlations, the 
possibility remains that the net repulsive interaction be- 
tween holes could be long-ranged. Such an unexpected 
result may indeed be connected with a true deconfinc- 
ment of spinon-like excitations for the kagome system, 
which would confirm speculation to this effect in 2dm& 
We comment further on the issue of deconfinement be- 
low. 



V. CONCLUSION 

We have illustrated a variety of novel phenomena asso- 
ciated with the effects of static, spinless impurities doped 
into the 5 = 1/2 antiferromagnetic Heisenberg model in 
kagome geometry. By analysis of spin-spin correlation 
functions and of the local magnetization patterns induced 
around impurities, we find that, despite the very short 
spin-spin correlation lengths of the highly frustrated sys- 
tem, dimer-dimer correlations develop over considerable 
distances due to a hole-induced freezing of their reso- 
nance. Furthermore, as a consequence of the continuum 
of states in the low-lying singlet manifold, there is no ev- 
idence for free local moments in the vicinity of impurity 
sites. 

Most importantly, we have shown from the total ener- 
gies of doped systems that impurities in the kagome lat- 
tice have an effective repulsive interaction, a phenomenon 
unique to our knowledge and in direct contradiction to 
the conventional expectation for short-range RVB sys- 
tems of a net pairing attraction. Furthermore, this effec- 
tive repulsion is found to be rather long-ranged in na- 
ture, a most unusual result which is presumably related 
with the very strong frustration inherent to the kagome 
geometry. 

Our comparisons of doping effects provide additional 
justification for the validity of the short-range dimer- 
RVB framework as a suitable description of the ground- 



state singlet manifold in the kagome antiferromagnet. We 
have shown that the dimer-RVB basis provides a valu- 
able means of probing the nature of spin and charge de- 
grees of freedom in a further class of quantum magnets. 
This realization of a highly resonant RVB state, arising 
as a consequence of the extremely frustrated nature of 
the system, also presents an important contrast to other 
low-dimensional spin systems which are known to be de- 
scribed by RVB-type wave functions, particularly with 
regard to the localization of free spins and to the pairing 
of holes. 

We summarize only briefly the consequences of our 
analysis for experiment. The majority of known kagome 
materials do not have complete filling of the lattice sites 
by a single value of the spin 5, and as such the study 
of impurities is a significant component of their charac- 
terization. Nonmagnetic impurities in the kagome lattice 
have the important property that they do not generate 
localized spins at the vacant site, and do generate by the 
nature of their extended range of influence a significant 
number of satellites for each peak. On heuristic grounds 
we believe that our analysis of the qualitative features to 
be expected in the 5=1/2 system are relevant also for 
5 = 3/2, and perhaps for higher half-integral spins. 

Finally, on the issue of deconfinement in highly frus- 
trated systems of dimensionality higher than one, we 
have shown clearly from the presence of the mutual re- 
pulsive interaction between holes that charge degrees of 
freedom are deconfined. Further, the absence of localized 
5=1/2 degrees of freedom around vacancies in the lat- 
tice demonstrates that spin and charge are deconfined. 
However, on the question of whether the spin excitations 
of the undoped system are appropriately described by 
deconfined spinon excitations, our results do not provide 
a proof of this hypothesis. 
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